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Many popular models of new physics beyond the Standard Model use a parity to stabilize weakly 
interacting, dark matter candidates. We examine the potential for the CERN Large Hadron Col- 
lider to distinguish models with parity stabilized dark matter from models in which the dark matter 
is stabilized by other symmetries. In this letter, we focus on signatures involving long-lived par- 
ticles and large amounts of missing transverse energy. To illustrate these signatures, we consider 
three models from the literature which are representative of a more general class of models with 
non-traditional stabilization symmetries. The most optimistic scenario can observe the proposed 
signature with a minimum of 10 fb _1 of integrated luminosity at design center of mass energy. It 
will probably take considerable longer to validate the stabilizing symmetry is not a simple parity. 
In all, we emphasize that the underlying symmetry that stabilizes weakly interacting dark matter 
has tremendous implications for the LHC and our understanding of the nature of dark matter. 
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There is an overwhelming amount of evidence for the 
existence of dark matter (DM). Numerous astrophysical, 
cosmological and direct detection experiments [l| provide 
a consensus picture: viable dark matter candidates must 
be stable, neutral under the Standard Model (SM), non- 
relativistic at redshifts of z ~ 3000 and generate the mea- 
sured relic abundance of h 2 fl DM = 0.1131 ± 0.0034 0. 
These requirements are general and upon implementa- 
tion provide a plethora of models with viable DM candi- 
dates. The most popular, e.g. supersymmetry (SUSY), 
Little Higgs and extra-dimensional scenarios, typically 
stabilize the dark matter with a generic parity (Z2) sym- 
metry 0, 0, H i, 0, H IH. Since particles are defined 
by how they transform under different symmetries, in 
essence only one type of DM candidate is being consid- 
ered! In this letter, we consider the observational conse- 
quences of dark matter stabilized with another symmetry 
(a "non-parity" symmetry) at the CERN Large Hadron 
Collider (LHC). We focus on a signal of two long-lived 
particles plus large amounts of missing transverse energy 
(fjr). 1 To demonstrate the viability of this signature, we 
devote the next section to show how such a signal is often 
suppressed for models with Z2 stabilized DM. We next 
briefly review three example models in which the DM is 
stabilized by a symmetry other than a Z2; there example 
processes that generate the signal are provided. Section 
III contains the specifics of our methodology including 
the SM backgrounds, detector and acceptance cuts. Also 
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1 In addition to long-lived particles, one can also study decay 
chains of particles charged under "non-parity" stabilization sym- 
metries. The result is specific kinematical distributions that are 
distinct from traditional parity stabilized models. We explore 
this research direction in the collaboration of [id . In a separate 
collaboration [Till , we focus on showing how the LHC can distin- 
guish strongly coupled hidden sectors with conformal or confined 
phases from perturbative ones using jets of missing energy. Such 
sectors can be the source of the dark matter candidates. 



included are the details of a fast detector simulation writ- 
ten to estimate the amount of missing energy generated 
from long-lived particles interacting with the ATLAS and 
CMS detectors. The relative suppression of Z2 stabilized 
models for the signal of two observed long-lived particles 
plus large amounts of missing transverse energy is a cen- 
tral point of this paper. This statement is quantified for 
an example effective theory in Section IV. The analysis 
therein includes the computing of the statistical signifi- 
cance for the "non-parity" signal processes from Section 
II. Also considered are the well-known SM backgrounds. 
Conclusions follow. 



I. SIGNAL SUPPRESSION FOR PARITY 
STABILIZED MODELS 

In this section we provide inductive reasoning to sup- 
port the assertion that models with Z2 stabilization sym- 
metries are either generally suppressed or distinguishable 
from models with non-traditional stabilization symme- 
tries. These arguments are valid for signatures of two 
long-lived particles plus f)r. To begin, consider the low 
energy effective theory of a model with parity stabilized 
dark matter, \- Below the cutoff A, we include the SM, 
X and the next-to-lightest parity odd particle, vb. We 
require vb to be charged under both the SM and Z2 sym- 
metries. We also require vb to have a lifetime long enough 
to ensure it will not decay inside of the detector yet is 
consistent with the strong bounds on charged, long-lived 
particles fl2l |. Now just by parity conservation, any pro- 
cess involving the production of parity odd particles must 
have an even number of such particles in the final state. 
In our scenario, this means a final state involving vb*vb, 
vb*ip vb*vb, ifi*4>X*Xi e tc. Now suppose a dedicated exper- 
imental analysis observes long-lived yb particles. From 
our effective description, we posit final state signatures 
with two observed yb particles plus frp are generally sup- 
pressed. Such a signature may be possible for models in 
which the dark matter is stabilized by other symmetries 
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such as Zi x Z 2 , Z 3 , SU(3) g \ ba\, etc. In Section III and 
IV, we will account for the relevant SM backgrounds. 

To see the signal suppression in the Z 2 case, note the 
most relevant operators involving couplings of tp with the 
SM are 



and 



Oif = gip c aAip 



2b = Mh ] htp*tp 



lb =gdtP*AiP 



(1) 



° 2f = x fctfc ^ (2) 



where A is the gluon, photon or Z boson, h is the higgs. 
Note, h is general and can represent extended higgs sec- 
tors such as two higgs doublets that are common in SUSY 
models. Oif (On,) are the operators involving fermionic 
(bosonic) tp. The same type of couplings are also possible 
for x with A the photon or Z boson. If both x an d tp are 
scalars, then the marginal operator 



3b = X 3 ip*ipx*X 



(3) 



is also possible. As well, if tp has the SM quantum num- 
bers of the higgs then 



04b = 



A 



O if = \ if rxh (4) 



is viable. If tp has SM charges which allow yukawa cou- 
plings to quarks and leptons, our assumption of a long- 
lived tp particle requires, As^, to be small for 



5 = X 5 tp c X l 



6 = A 6 ip c XQ- 



There may be additional higher-dimensional operators 
that may permit decay of tp into the SM. We assume A 
is sufficiently large to suppress the effects of these oper- 
ators. 

One way to generate the signal with our effective de- 
scription is through pair production of the tp particles. 
One of the ip particles can subsequently emit two x par- 
ticles via equation [3l It is also possible for tp to emit 
a SM higgs boson or an off-shell Z or photon via equa- 
tions [1] and [2 These would then need to decay into two x 
particles. The first (second) scenario is suppressed by the 
amount tp (as well as Z*/j*) goes off-shell. The diagrams 
involving a virtual SM higgs, ho, have additional off-shell 
suppression if m/ l0 < 2m x . To get a numerical estimate 
of the effect of off-shell suppression, note the integral over 
a Breit-Wigner resonance for one virtual particle, V, is 



T V M V 



(6) 



where 8 



max(min) 



= tan 



-Him 2 



max(min) 



M 2 )/T V M V ). 

Here m max and m m i n are the appropriate limits of the 
phase space integration (l3j . In the limit where the vir- 
tual particle is on shell max — #min - * T. As an example, 
for the process involving equation [31 the integral goes as 
#max - 6> min ~ 0(0.01) for a 500 GeV ip and 100 GeV X - 



The SM higgses in the diagrams described above can have 
heavy masses, 1 TeV > rrih > 2m x , which may slightly 
reduce the amount of off-shell suppression; in this case 
often mh > 2 mz and there is an additional suppression 
from the SM higgses preferentially decaying into W and 
Z bosons. Note the partial width for Hq — » XX ls 



* he — 'XX scalar 
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16nmh<_ 

A K 2 HT-ha 

32^ 



1 — Am 2 /r 



1 - Am 2 x /m 2 ho 



(7) 



3/2 



(8) 



where K\ (k 2 ) is the coupling for scalar (fermionic) X- 
Here A = 2, 4 for Weyl and Dirac fermions, respectively. 
The partial width for decays into W and Z bosons goes 
as r ~ GFrn\ where Gf is Fermi's constant. Thus, as 
TOft increases, the higgs decays more frequently into W 
and Z pairs due to the m\ enhancement. For example, a 
1 TeV (300 GeV) higgs decays into 100 GeV scalar (Dirac 
fermion) DM 6 x 10 -2 % (11%) of the time. Here we have 
assumed Ki = 0.5. Even when the SM higgs decays at 
the 10% level, an experimentalist can search for ip*ijj + ho 
where ho decays to the SM. As we will discuss in Section 
III, tagging long-lived particles is relatively easy. Thus, 
precise measurement of the SM higgs branching fractions 
can uncover an invisible decay width. An experimentalist 
can potentially identify if a pr signature comes from SM 
higgs decay. 

Another way to generate the ip*ip + fir signature is 
through pair producing charged higgses. Such higgses 
are common in SUSY models. The higgses decay on- 

f J7ty. 



(5) shell via the operators in equation [4] for rrih > 



Now, in two higgs doublet models, the partial width for 
each higgs decay channel into the SM goes as T ~ Gf m\. 
Thus, from the arguments in the previous paragraph, the 
h — > ipx* decay may be suppressed. Even with a signif- 
icant partial width, the experimentalist can still largely 
determine if higgs decay is responsible for the tp*ip + fir 
signature. The higgs decay width into tp x* can be mea- 
sured by comparing pp — * h* h — * tp + SM + fi)r events to 
the signal. In particular, this process allows the recon- 
struction of both of the virtual higgses. The critcal point: 
If an experiment observes the tp + SM + fir signature and 
the SM daughters can be reconstructed to the higgs, it 
is observationally clear that higgs is parity even. If the 
mass of the other virtual higgs can be reconstructed, it 
is clear we have the case where parity even states can 
produce the ip*tp + fir signature. The higgs mass from 
the h — > ipx* decay can be reconstructed with a trans- 
verse invariant mass. This is in analogy to W — > I v mass 
reconstruction in WW decays [l4j]. As before, observing 
the higgs decay into the SM daughters is not as problem- 
atic because the tagged tp particles greatly reduce the SM 
background. The transverse mass is defined as 

M T = V(E T + fir) 2 -(PT+0T) 2 (9) 

where E T = pS^+M^ and pr is the transverse momentum 
of the tp particle. In addition, we take jfr — ^2pt visible 
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FIG. 1: An example reconstructed transverse mass for pp — » 
h* h —tip x* + 1 9 with 500 GeV charged higgs decaying into 
a 200 GeV long-lived tp and 100 GeV DM y*. Th e most 
important aspect of this plot is the kinematic edge at 500 
GeV. We apply the cuts consistent with ATLAS and CMS 
collaborations with a transverse momentum cut of 200 GeV 
and the cuts of equations 1161 and 1191 



and ^3p = Pt 2 ■ As an example, we simulate the process 
pp — > h* h — > tpx* + SM for a 500 GeV charged higgs 
and plot the reconstructed transverse invariant mass in 
Figure 1. Note the transverse mass has a kinematic edge 
at the given higgs' mass. 

The careful reader may argue that adding additional 
particles to the effective theory will generate large miss- 
ing energy signatures for models with Zi stabilized dark 
matter. There are three types of new particles which can 
produce non-trivial consequences: a particle, (p, charged 
solely under the SM, the next heaviest parity odd par- 
ticle, tp', which is also charged under the SM and <p s 
a singlet under the SM and the Z2 symmetry. In the 
first case, the coupling A7 (p* tpx 1S allowed; thus, gener- 
ating the tp*ip + fjT signature with the minimal amount 
of suppression mirrors the discussion for the production 
and decay of charged higgses via equation |4j We note <f> 
cannot have random SM quantum numbers which could 
forbid yukawa couplings with SM particles. This is gener- 
ally true because (1) (p has the quantum numbers of ip°x* 
and (2) tp must decay through relevant or irrelevant oper- 
ators to y an d the SM to avoid bounds on stable charged 
particles [lj]. Like the discussion of the two higgs dou- 
blets, <p should decay more prominently to the lighter SM 
particles than tp and y; therefore, the tp*tp + frp signal 
may be suppressed. As before, if there is a non-trivial 
decay width into tp y* , plotting the transverse invariant 
mass could provide identification of tp*tjj + pr production 
from a parity even (p. Consider now two new effective 
descriptions where we separately add tp' and <p s to the 
original effective theory of the SM, tp and y. In the tp' 
effective theory, it is clear that there are no tp'tpx cou- 
plings; however, if tp and tp' have the same charges under 
the SM, they can mix by an angle, 9, that is a function 
of the mass ratio m$ / m^i . Any cross section involving 



pair production of a heavy tp' with a tp* tp y* y final state 
is suppressed by at least 9 4 . We have not demanded 
that tp' is long lived and it can conceivably decay to the 
SM and y via operators analogous to equation [H thus, a 
small branching fraction into our signal may be possible. 
In the <p s effective theory, <p s must have 0(1) couplings 
to the SM in order to be produced at the LHC; thus, it 
generally mirrors the SM higgs and charged higgs discus- 
sions above. The signal could be generated from <p s pair 
production and decay via (p s tp*tp and <j> 8 y*y couplings. 
Like before, the parity charge of <p s can be ascertained 
by considering the <p s (p s — > y* y + SM process. Both cp s 
masses could be reconstructed in analogy to ZZ decay 
to two neutrinos and two visible SM particles [H^ . The 
reader might now argue that a new symmetry or large 
coupling could be implemented so that <ps decays solely 
to tp and DM pairs. Likewise, a similar argument could 
be constructed so <p decays solely to DM and tp. The 
large coupling case implies a new accidental symmetry. 
In both cases, tp and y would then have to be charged 
under these new symmetries. This is tantamount to the 
new stabilization symmetry we sought to distinguish in 
the first place. 

An effective theory with some number of <p, tp' and 
4> s particles is also possible. Given all of the arguments 
above, we assert such an effective theory either gener- 
ates a suppressed signal or can be identified as a sce- 
nario with a Z2 stabilizing symmetry. If further convinc- 
ing is necessary, one can list all of the tree-level decay 
topologies which result in a tp* tp y* y final state where 
the mother particles have non-specific charges. The least 
suppressed diagrams from this exercise are analogous to 
examples above. As a final note, since both <p and <p s 
can potentially have tree-level SM couplings, generally 
precision electroweak measurements place constraints on 
their mass to roughly 5 TeV [13] ■ In the next section, 
we briefly outline three models in the literature in which 
the dark matter is stabilized by a non-traditional "non- 
parity" symmetry. We will subsequently show how these 
models generate relatively large tp*tp + f}r signatures. Ul- 
timately, our signal will be a statistically significant ex- 
cess of missing energy in this channel. We present each 
model's results in comparison to an effective theory with 
a Z2 stabilization symmetry This comparison effective 
theory has a particle content of the SM, the SM higgs, tp, 
tp' and Dirac fermion DM. The missing energy is gener- 
ated from the process h — > x*X f° r a SM higgs mass of 
350 GeV. The DM is considered to be Dirac DM in order 
to maximally enhance the invisible branching fraction. 



II. EXAMPLE MODELS 

In this section, we briefly outline the relevant aspects of 
three models in the literature that can generate our sig- 
nature. The first model stabilizes the dark matter with a 
linearly realized discrete symmetry generated from spon- 
taneous symmetry breaking. The second stabilizes with 
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a i?3 symmetry generated from baryon number. The fi- 
nal model has DM stabilized by a parity symmetry. A 
light, dark gauge boson which couples to the DM helps 
to mediate DM annihilation. 



A. A Companion Model 

In a companion paper, we present a class of models in 
which spontaneous symmetry breaking is used as a mech- 
anism to naturally correlate the weak scale and the mass 
scale associated with the dark matter annihilation cross 
section. See for details. The featured model has a 
SU(2)d that commutes with the SM; the symmetry is 
"broken" along with electroweak symmetry through an 
extended higgs sector to generate a linearly realized Zi 
symmetry. The result is new dark gauge bosons with a 
mass of order the weak scale. These models can accom- 
modate additional messenger vector- like heavy fermions, 
Q, whose mass is, for our purposes, a relatively free pa- 
rameter [l2j]. Q is charged under the SM as well as the 
SU(2)d- For our signature, we consider these vector- like 
fermions, Q, as long long lived. We are interested in the 
effective operators 



T =g'Q*aAQ 



(10) 



where A is the new gauge bosons that does not couple di- 
rectly to the SM. As shown in [ijj], the new gauge bosons 
decay into two DM candidates, Xi with a 100% branching 
fraction via the operator 



08 = g'x'oAx 



(11) 



Thus, the signal for this model at lowest order is pp — ► 
QQ*A — > QQ*pr- There is a seemingly irreducible 
background of pp — > QQ*. We will show how to dis- 
tinguish this background in the following sections. Fur- 
ther, we also show a large excess of missing energy in this 
channel in comparison to the Z2 scenario described in the 
previous section. In the coming analysis, we refer to Q 
and A as ip and x, respectively, in order to be consistent 
with the notation in previous sections. 



B. Agashe and Servant 

Agashe and Servant (AS) proposed a class of grand 
unification models within the framework of warped extra 
dimensional space-time. Please see [l8[ for more details. 
In these models, the SM along with new exotic matter 
fills out the 5*0(10) spinor representations. Each spinor 
multiplet is charged under baryon number. Notably, all 
of the non-SM particles transform under a -Z3 symmetry 
generated from the particles' color and fractional baryon 
number 



rj —>■ exp 



2ni B 



(12) 



Here rj is a generic of any particle in the theory and £?, n c 
(n c ) are the particle's baryon number and number of col- 
ors (anti-colors). A natural DM candidate, vr, with the 
SM charges of a right-handed neutrino, is in same the 
spinor multiplet as the SM right-handed top, tR, with 
B = —1/3. Also of interest to us is an exotic, SU(2)l 
doublet fermion, Ql, which transforms as the fundamen- 
tal under color SU(3). Ql is charged under the Z 3 and 
is in the same spinor multiplet as vr; thus, it can be 
relatively light compared to the other new exotic states. 
Because the 50(10) is broken to the SM gauge group, 
there are additional gauge bosons beyond the SM that 
are charged under the SM and Z3 symmetries. The X 
boson has the same SM quantum numbers as Ql but with 
a different baryon number of B = —2/3. We note because 
the SO(10) was broken with varying boundary conditions 
than the SM, proton decay constraints are eliminated. 
We focus on the coupling 



Og = g^Q* L av R X 



(13) 



where g§ is the Pati-Salam gauge coupling. There is an 
additional Pati-Salam gauge boson, X s , which is impor- 
tant when considering the decay of Ql to the SM and 
vr. This Ql decay can only occur through equation 1131 
where the X boson mixes with the X s boson. The X s 
boson subsequently decays to vr and tR via 



O, 



g 5 t* R a vr X s 



(14) 



The mixing between the X and X s bosons occurs through 
a vaccum expectation value (vev) that can be effectively 
dialed. We assume a value that generates long-lived 
particles consistent with pjj]. Consider the AS model 
where the SO(10) gauge group is broken on the infrared 
brane. In that scenario, the X bosons can be relatively 
light. 2 In the case where the SO(10) is broken on the 
UV brane, the X bosons must be heavier than about 
3 TeV. We consider the former and take the process 
pp — ► Ql Q* l ~ * vrv* r X X* with X being long lived as 
our signal. We consider only QCD production involving 



0n = S s QIoGQl 



(15) 



where G is the gluon. We take the mass hierarchy of 
mQ L > mx + m VR . The signal is enhanced because Q l 
decays to X* vr with a 100% branching fraction. In the 
analysis section, we refer to X and vr as ip and Xi re ~ 
spectively, in order to be consistent with the notation in 
previous sections. 



C. Light Hidden Dark Sectors 

Pospelov, Ritz and Voloshin proposed a class of mod- 
els in which the DM has "secluded" interactions with 



2 We thank K. Agashe for alerting us to this point. 
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the SM. See [2l[ for details. This is done so their mod- 
els can be consistent with increasingly stringent bounds 
from elastic DM scattering with nuclei in direct (as well 
as some indirect) detection experiments. The DM anni- 
hilates through metastable mediators which in turn even- 
tually decay to the SM. All of the secluded versions of 
their models have the characteristic that the DM mass 
is larger than the mediator mass. As an example, they 
construct a model in which a new U(l) gauge boson ki- 
netically mixes with SM hypercharge. After spontaneous 
symmetry breaking, the new gauge boson gets a mass and 
has a lifetime which, for most of the parameter space in 
the model, is stable on detector time scales. 3 This model 
differs from the model in Subsection IIA in that the for- 
mer uses spontaneous symmetry breaking to generate the 
dark matter stabilization symmetry. 

More recently, Arkani-Hamed, Finkbeiner, Slatyer and 
Weiner (AFSW) proposed a class of models to provide 
a dark matter interpretation for the PAMELA anomaly. 
The PAMELA collaboration searches for indirect signals 
of dark matter by observing high-energy cosmic rays from 
orbit. The collaboration observes an anomalously large 
number of positrons in the cosmic ray fluxes which do 
not conform to the theoretically expected background. 
In constrast, the anti-proton to proton flux ratio does 
conform to theoretically expected background 20]. The 
authors proposed a class of models where a new force 
helps to mediate the DM annihilation. See [22j for de- 
tails. To be consistent with the positron and anti-proton 
measurements, the new gauge boson(s) must have a mass 
less than twice the mass of the proton. Thus, kincmat- 
ically, the new gauge bosons therefore can only decay 
into the lighter leptons. The fact that this new gauge 
boson couples to leptons illustrates a primary difference 
between this and the companion model in the first sub- 
section. 

For both scenarios, we assume a new light gauge boson 
with a mass of 1 GeV. The gauge boson does not decay 
inside the detector and constitutes the missing energy. 
This is in analogy to the signal process in Subsection 
IIA. For the AFSW model, this means we assume the 
new gauge boson to be long lived and not decay into 
leptons inside the detector. 4 We also assume a vector-like 
quark production mechanism via the effective operator in 
equation [TUl 



III. METHODOLOGY 

A huge key to our signature are the particles, ip, which 
are long lived and charged under both the SM and dark 
matter stabilization symmetry (-ies). Here we briefly re- 



3 We thank A. Nelson for this observation. 

4 We thank N. Weiner for suggesting this point for the model of 
Arkani-Hamed, et al. 



view some important aspects of the collider physics of 
long-lived particles. A review can be found in [231 ]. 

If V' is charged under SU(3) color, it, after production, 
forms color singlet states with valence SM quarks. These 
■0-hadrons interact with the detector by exchanging SM 
quark partons with the material therein. This nuclear 
scattering process causes these ^-hadrons to potentially 
charge flip from charged to neutral (and vice- versa), lose 
energy and/or even stop inside the detector. The heavy 
tj) parton inside the V'-hadron carries a dominant frac- 
tion of the total energy of the system with the interact- 
ing SM quark partons carrying approximately O(0. 001) 
fraction of the rest. Thus, the nuclear interactions of 
these ■f/'-hadrons are dominantly low-energy QCD with 
the changes in momentum due to interactions with the 
detector small relative to the total momentum. If tp is not 
colored, then it interacts only electromagnetically with 
the detector. In this and in the charged ?/>-hadron case, a 
highly ionized track is produced in the detector's tracking 
chamber. Upon reaching the muon detector the charged 
particles look like a slow, heavy "muon." 

We are interested in large frr signatures generated 
by the dark matter candidates. Since the hadrons can 
potentially charge flip to neutral as well as lose energy 
through interactions with the detector, it is important to 
quantify the amount of missing energy in each event. To 
address this, we wrote a fast detector simulation which 
parametrizes the charge flipping and energy loss effects 
accounted for in GEANT3 and 4 (24|. In our simple 
setup, we assume perfect lead calorimeters. The hadron 
interactions therein are assumed to be perfectly opaque 
while the signals of the interactions are perfectly translu- 
cent. The detector response is extrapolated from SUSY 
gluino R- hadron nuclear interactions for iron (2f| . We do 
not include ionization effects in our simulation; such dom- 
inates over nuclear interactions only when the hadrons 
are at very small velocities. The geometry for the AT- 
LAS and CMS detectors in terms of nuclear interaction 
lengths is taken from [261 ]. In Figure 2a, we plot the 
missing energy generated from two tagged tp particles af- 
ter interacting with the calorimeters in our ATLAS and 
CMS detector simulation. The ip particles have mass of 
300 GeV and transform as a fundamental under color 
SU(3). The similarity between the two curves is due 
to the similarity of the ATLAS and CMS calorimeters 
and inner detectors in nuclear interaction lengths. It is 
comforting that the resulting distribution favorably ap- 
proximates the ATLAS GEANT3 simulation for stable 
gluinos in Figure 9 of 27]. In Figure 2b, we plot the 
missing energy for the pp — > + fr signal as well as 
the pp — > tp*tp background for 10 fb" 1 of luminosity. We 
used the companion model scenario described in the pre- 
vious section as well as the parameters listed in Section 
IV. We also used the ATLAS and CMS cuts described 
in the next paragraph. Here the difference between the 
ATLAS and CMS curves is due to the differing kinematic 
cuts. It is clear the missing energy generated from the 
long-lived particles traveling through the detector has the 
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FIG. 2: (a) Events vs. missing energy {frr) for two 300 GeV 
quarks after transversing the simulated ATLAS (dashed) and 
CMS (solid) inner detector and calorimeters, (b) The same 
plot as (a) but with the kinematic cuts from equations ll6ll9l 
Also shown for ATLAS (dot-dashed) and CMS (dotted) is the 
tp'ip + f)r signal in addition to the ip*ip background for the 
companion model in section IIA. g' is set to be electroweak 
strength, (c) The normalized cross section versus velocity for 
a 300 (solid), 900 (dashed) and 2000 (dot-dashed) GeV long- 
lived vector-like quarks in the central detector region. 



potential to obscure the signal. This plot demonstrates 
the necessity of writing the fast detector simulation. 

Any SM process which can generate muon tracks 
and/or hard jets that punch through to the muon de- 
tector provides a potential background to our signal. To 
begin, it should be noted the (near-)massless particles 



have a velocity, (3, near speed of light (j3 = 1); the heav- 
ier hadrons travel at much slower velocities and arrive 
at the muon detector "out of time" with the rest of the 
event. In figure 2c, we plot the velocity for hadrons of 
different masses. It is clear most of the events have j3 < 1. 
An additional complication occurs if the long-lived par- 
ticle arrives too late. Each proton beam at the LHC is 
collimated in bunches and is scheduled to be nominally 
separated by 25 nanoseconds [28[. Very slow hadrons 
could be mis-reconstructed with the wrong event. Fur- 
ther, without the proper cut, cosmic ray muons could 
provide an additional background. Because up to three 
events could be in the detector at one time, the larger 
ATLAS detector requires a higher lower velocity cut. We 
adopt velocity cuts suggested by the ATLAS and 
CMS collaborations 

0.7 < /?atlas < 0.9 and 0.5 < (3 C ms < 0.8. (16) 

To limit the effect of hard jets punching through to the 
muon detector and faking the long-lived particles, the 
ATLAS collaboration adopted additional stringent pr 

cuts m 

Pt atlas > 250 GeV and p T C ms > 30 GeV. (17) 

For the ATLAS (CMS) cut, at least one (both) of the 
"muons" must satisfy the cut. The smaller CMS px cut 
is due to the collaboration requiring both a highly ioniz- 
ing track in the inner detector as well as muon detector 
information to eliminate SM backgrounds. This ionized 
track is subject to the velocity cut in equation[T|)J Given 
the measured velocity and momentum, CMS also requires 
the average reconstructed particle's mass, m, to be 



m C MS > 100 GeV. 



(18) 



where P — (3m/ yl — j3 2 . Full simulation of ionization 
effects are beyond the scope of this work. As a simplifica- 
tion, we simply retained all events with charged hadrons 
in the simulated CMS tracker. For the simulated events 
we tag at least two long-lived particles in the muon detec- 
tor. 5 We require these particles to be in the triggerable 
part of the muon detector 



M < 2-4. 



(19) 



ATLAS requires no hard jets in the calorimeters to come 
within a cone of Ai?ATLAS < 0.4 of the "muon" track. 
CMS requires the track from the inner detector to be in 
a Ai?civis < 0.1 cone of the "muon" track in the muon 
detector. Here AR = \J Arf + A(j) 2 where rj (4>) is the 
pseudorapidity (transverse angle) of the event. For the 



5 If an odd number of long-lived, visible particles are detected in 
addition to the i/>*i/> + frp signature, it is clear a symmetry other 
than a parity is responsible. 
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ATLAS cuts, the collaboration finds a signal to back- 
ground ratio of S/B = 2.6 x 10 3 for an integrated lumi- 
nosity of 1 fb _1 . Similarly, for the CMS cuts, that collab- 
oration find a "background free region" for 100 pb _1 [29| . 
With such a high expected S/ B ratio from both collabo- 
rations, our analysis will not focus on the traditional SM 
backgrounds. 

Because we are looking for signatures with at least 
two long-lived particles plus large ftr, there are irre- 
ducible backgrounds to consider. To begin, the Z bo- 
son can decay to neutrinos and generate the signal via 
pp — ► xjj*ip Z — ► ip* ipv* v. The precise Z invisible width 
is well known and can be accounted for. Further ver- 
ification the Z is the cause of the missing energy can 
come from counting the number of pp — > ^*ip Z — > ip*ip 1 1 
or — > ip*ipqq events and reconstructing the Z invariant 
mass for identification. To be conservative, in addition 
to the invisible Z width, we add an invisible SM higgs 
decay width. We assume the operator, k h x*X? where 
k ~ 0.5. x is taken to be Dirac dark matter to maximize 
the invisible decay width. See equation [51 We require 
a statistically significant number of signal events above 
these irreducible backgrounds. Finally, especially for the 
first and third models in the previous section, pair pro- 
duction for ip particles dominates typically by orders of 
magnitude over the ifi*ip + frp signature. The tp* ij) pro- 
duction is generally back-to-back. To further define the 
signal, we require the angle in the transverse plane to be 

COS0 > -0.9. (20) 

to effectively eliminate the ip*ip production. 

IV. ANALYSIS AND RESULTS 

In this section, we show a statistically significant excess 
of missing energy for our signal of two long-lived parti- 
cles, tp, plus large amounts of missing transverse energy, 
pr- To characterize this signal, we focus on the three 
models described in Section II. All events are simulated 
for the LHC at design 14 TeV center of mass energy. 
Our simulations use CTEQ Set 4M parton distribution 
functions. a s is calculated at two loops with the renor- 
malization and factorization scales set by vi/2. After 
implementing the cuts described in the previous section, 
we assume 100% tagging efficiency for the visible heavy, 
long-lived particles. The masses of the stable particles 
are taken as 

my = 300, 600 GeV. (21) 

The mass of the DM candidates is set to be 

m x = 100 GeV. (22) 

In the first section, we argued models with parity sta- 
bilized dark matter generated suppressed + fir sig- 
natures. To quantify this statement we consider a par- 
ity conserving effective theory. As discussed before, we 



Background Events 



ATLAS 


300 GeV V|600 GeV V 


300 GeV V 




(with cos <f> cut) 


(without cos <j> cut) 


Z —* u*u 


71.2 


6.1 


130 


ho -* X*X 


0.126 


0.01 


0.183 


CMS 








Z — » v*v 


49 


5.4 


73.6 


ho -^X*X 


0.055 


0.006 


0.066 



TABLE I: The simulated background for 10 fb" 1 of integrated 
luminosity. Here the Z —* v*v and ho — > X*X labels corre- 
spond to the background processes: pp — » ip*^> Z — -> ip*^jv*v 
and pp — > ip*t[>h — * ip*ipX*X- The former is the well known 
invisible Z decay; the latter is the assumed invisible higgs de- 
cay width for the parity odd effective theory described below. 
X is taken to be a 100 GeV Dirac DM. The larger suppression 
ho — * XX i s due primarily to the heavier 300 GeV SM higgs 
forcing the i[> off-shell as well as the 11% invisible branching 
fraction. 

choose the effective theory to have the SM, SM higgs and 
parity odd ip and Dirac DM, Xj as the relevant low energy 
degrees of freedom below the cutoff A. We require the 
SM higgs to have a significant invisible branching frac- 
tion into the DM. Conservatively, we take the coupling 
from equation [2] to be A2/ = 2 with A = 1 TeV. As dis- 
cussed below equation [SJ a 300 GeV SM higgs decaying 
into Dirac DM has a relatively enhanced invisible width. 
For these parameters, the generated background events 
are listed in Table I. As shown, the process involving the 
invisible higgs decay width is generally down by a fac- 
tor of 100 from the process involving a Z invisible width. 
This is primarily due to the amount the long-lived tp must 
go off-shell to accommodate the decay. Further, since the 
SM higgs couples to all of the massive SM particles and 
must be heavier than twice the DM mass in order for 
the process to avoid additional off-shell suppression, we 
remind the reader this suppression is a general feature. 
We now focus solely on the invisible Z background. 



A. A Companion Model 

The dark gauge coupling in the companion model is 
taken to have the same strength as the electroweak cou- 
pling measured at the Z mass g' = g ew (Mz)- In Figure 
3a, we plot the cosine of the angle between the recon- 
structed i/)-hadron momentum in the plane transverse to 
the beam direction. This is done for the pp — > ip*%p 
background as well as the pp — > ip*ipX signal for both 
the ATLAS and CMS detectors. We remind the reader 
X (tp) is the new heavy gauge boson (vector-like quarks) 
described in Section IIA. It is clear the missing energy 
signal is swamped and most of the events are back-to- 
back in the transverse plane. Please also refer to Figure 
2b; there we plotted the reconstructed missing energy 
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O 100 200 300 400 500 600 

Missing Energy (GeV) 

FIG. 3: (a) Events vs. cos<£ for ATLAS (solid) and CMS 
(dashed) for the %p*ij) background along with the ip*ip fir sig- 
nal. The missing energy is from the emission of a 100 GeV 
heavy gauge boson as described in Section IIA. The tp mass 
is 300 GeV. (b) The same plot as the first but rendered to 
illustrate the signal events, (c) For the model in Section IIA, 
signal events vs. f)r with the cos</> < —0.9 cut for ATLAS 
(solid) and CMS (dotted). All the plots implement the kine- 
matic cuts in Section III. 



from the ^-hadrons. The missing energy signature is 
also swamped. In Figure 3b, we redisplay the same plot 
to emphasize the signal events. The dark matter "kicks" 
the reconstructed ?/>-hadrons from being exactly back-to- 
back. After applying the coscj) cut, equation l20l none 
of the pp — » events remained. We plot the resid- 
ual signal missing energy events in Figure 3c for ATLAS 



Companion Model 

ATLAS CMS 







S/B 


s/Vb 


S/B 


300 GeV ip 


34.3 


4.1 


29.4 


3.92 


600 GeV ip 


10.1 


4.1 


9.5 


4.1 



Agashe- Servant 



ATLAS 



CMS 





s/Vb 


S/B 


s/Vb 


S/B 


300 GeV V 


34.3 


4.1 


29.4 


3.92 



Light Hidden Sectors 



ATLAS 



CMS 





s/Vb 


S/B 


S/VB 


S/B 


300 GeV %p 
600 GeV ip 


5.33 
0.75 


0.63 
0.30 


8.76 
1.08 


1.04 
0.47 



TABLE II: The signal significance and signal-to-background 
ratio for the described models in Section II for 10 fb _1 of 
integrated luminosity. Included is the dominant invisible Z 
background. 



(solid) and CMS (dotted). The generated signal events 
for 300 GeV long-lived ip and 100 GeV DM for 10 fb" 1 
of integrated luminosity is 



Satlas = 290 



Scms = 192. 



(23) 



For a 600 GeV long-lived ip with the same luminosity, 
the generated signal cross section is 



S'a 



TLAS 



25 



S, 



CMS 



22. 



(24) 



In Table II, we summarize the signal significance as well 
as the signal-to-background ratio for all three example 
models. For the companion model, it is clear the signal 
can be observed for the chosen parameters. 



B. Agashe and Servant 

In this analysis, we take the Pati-Salam coupling to be 
55 = 0.7 with rriQ L — 600 GeV. The long-lived particle, 
if), is the heavy Pati-Salam gauge boson with a 300 GeV 
mass; recall, this X boson is a color triplet and picks up 
valence quarks to form the long-lived ^-hadron. Because 
Ql decays immediately, in this scenario there is no pp — > 
QlQl background with which to compete; also note, 
the Pati-Salam gauge boson-gluon coupling is through 
a suppressed dimension six operator. Thus, there is no 
competition from pp — > ip*ip as well. We still use the 
cos4> cut to reduce the Z — > v*v background. Applying 
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g 
S 




100 200 300 400 500 
Missing Energy (GeV) 




100 200 300 400 500 600 

Missing Energy (GeV) 

FIG. 4: Signal events vs. $r for ATLAS (solid) and CMS 
(dotted) in the AS scenario described in Section IIB (a) and 
the light hidden dark scenario described in Section IIC. (b) 
In both plots, the cos<j) < —0.9 cut is taken in addition to all 
of the kinematic cuts in Section III. 



all of the cuts listed in the previous section, we find signal 
events for 10 fb _1 of integrated luminosity of 



Satlas = 12 



S, 



CMS 



= 38. 



(25) 



In Figure 4a, we plot the missing energy for the signal 
for both ATLAS (solid) and CMS (dotted). In Table II, 
the signal significance and signal to background ratio is 
calculated. 



much as in the companion model. This increases the 
overall cross-section. The larger cross section is down by 
a factor of a — </ 2 /47r and is not competitive with the 
pp — > background. Although a 1 GeV light gauge 
boson does not generally "kick" the long-lived particles 
as much as the heavier dark gauge boson from the model 
in Section IIA, we still found it better to employ the cos (j> 
cut in order to reduce some of the SM background. For 
300 GeV ip and 10 fb _1 the generated signal events are 



Satlas = 45 
For 600 GeV ip 

Satlas = 1-85 



Scms = 74. 



Scms = 2.52. 



(26) 



(27) 



In Figure 4b, we plot the resulting missing energy. The 
signal significance and the signal-to-background ratio is 
summarized for ATLAS and CMS in Table II. A statisti- 
cally significant signal should be possible for 100 fb _1 of 
integrated luminosity for the case when — 300 GeV. 



V. SUMMARY 

In this letter we argued the signal of two long-lived 
particles plus large missing transverse energy were sup- 
pressed for models in which the dark matter and the long- 
lived particles were stabilized by a common parity (Z2) 
symmetry. Models stabilized by a "non-parity" symme- 
try may not be suppressed for this signal. We used this 
process to distinguish between the two scenarios. Specif- 
ically, we analyzed this signal process for three models in 
the literature which are stabilized by symmetries other 
than a parity. For each we showed statistically significant 
signals for our choices in parameter space. By compari- 
son, an example parity symmetric effective theory gener- 
ally contributed negligibly. Finally, to generate reason- 
able estimates of the amount of missing energy generated 
by the missing dark matter candidates, we wrote a fast 
detector simulation to account for the energy loss effect 
for heavy long-lived particles interacting in the ATLAS 
and CMS detectors. 



C. Light Hidden Dark Sectors 

For this model, we repeat the same analysis as in the 
first subsection but with a 1 GeV long-lived gauge boson. 
Again, we have assumed the hidden gauge coupling to 
be electroweak. Because the dark gauge boson is much 
lighter, it does not force the i\) particle go off-shell as 
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